IC-SAD-90/35E-1 945  Effects  of  Frequency,  Temporal,  and  Spatial 

Averaging  on  Image  Interference 


David  Taylor  Research  Center 


Bethesda,  Maryland  20084-5000 


AD-A222  501 


DTRC-S  AD-90/35E-1 945  February  1990 

Ship  Acoustics  Department 
Departmental  Report 

Effects  of  Frequency,  Temporal, 
and  Spatial  Averaging  on  Image 
Interference 


Paul  T.  Arveson 


DT1C 

ELECT  E  i 
MAY241990 


I 


Approval  for  public  release; 
distribution  is  unlimited. 


90  05  23  04  5 


» 


/ 


MAJOR  DTRC  ^BCHNICAL  COMPONENTS 


CODE  01 1  DIRECTOR  OF  TECHNOLOGY,  PLANS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


This  document  contains  information  affecting  the  national  defense  of  the  United  States 
within  the  meaning  of  the  Espionage  Laws,  Title  18,  U.S.C.,  Sections  793  and  794.  The 
transmission  or  revelation  of  its  contents  in  any  manner  to  an  unauthorized  person  is 
prohibited  by  law. 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1 .  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semHormal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case-by-case  basis. 


NOW-OTNSROC  5602/50  (R«v.  2-86) 


CONTENTS 


Page 


ABSTRACT  .  1 

ADMINISTRATIVE  INFORMATION  .  1 

ACKNOWLEDGEMENTS  .  1 

INTRODUCTION  .  2 

AVERAGE  OVER  A  FREQUENCY  BAND . 3 

EFFECT  OF  ROUGH  SEA  SURFACE . 4 

AVERAGE  OVER  A  FINITE-SIZE  SOURCE . 8 

APPENDIX  A.  DERIVATION  OF  IMAGE  INTERFERENCE 

TERM  FOR  A  BAND  OF  FREQUENCIES . 11 

APPENDIX  B.  COMPUTER  PROGRAM  LISTING.  .  . . 13 

REFERENCES . 19 


FIGURES 


1.  Geometry  of  the  Image  interference  effect.  ...  2 

2.  Detector  output  level  versus  geometry . 6 


j  Aooession  For  ,  j 

NTIS  GRAJel 
DTIC  TAB 
Unannounced 
Justif lcatic 

% 

□ 

n 

Ry 

Distribution/ 

Availability  Codes 

Dist 

ti 

Avail 

Spec 

and/or 

lal 

ABSTRACT 


This  report  describes  a  significant  phenomenon 
affecting  the  propagation  of  underwater  sound  at 
ranges  less  than  1000  yards:  the  image  interference 
or  Lloyd-mirror  effect.  Previous  models  for  the 
effect  do  not  adequately  take  into  account  the 
effects  of  frequency  averaging  and  spatial  averaging 
due  to  sea  surface  roughness  and  finite  source  size. 
These  refinements  on  the  basic  theory  are  developed 
here.  Also  a  clarification  of  the  meaning  of  the 
surface  reflection  coefficient  is  given  for  the  case 
of  omnidirectional  sources  and  receivers  —  the 
situation  that  prevails  in  ship  radiated-noise 
trials . 
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David  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC) .  The 
funding  came  from  several  projects  sponsored  by  the  Naval  Sea  Systems 
Command  (NAVSEA  55N) . 
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INTRODUCTION 


The  reflection  and  scattering  of  underwater  sound  from  the  ocean  surface 
is  a  phenomenon  that  has  been  studied  extensively  since  the  time  of  Rayleigh 
(1).  In  fact,  the  approach  taken  by  Rayleigh  has  been  followed  by  many  invest¬ 
igators  up  to  the  present  time.  Most  of  the  theoretical  literature  deals  with 
the  Rayleigh  problem:  a  plane  wave  insonifying  a  small  area  on  the  sea  surface. 
Reports  on  experiments  often  deal  with  an  approximation  to  this;  that  is,  a 
narrow-beam  transmitter  and  receiver  aimed  at  a  small  area  on  the  sea  surface 
or  some  model  of  the  surface. 

The  situation  of  interest  in  ship  acoustical  trials  is  best  represented  by 
an  omnidirectional  source  and  receiver  under  a  moderately  rough  sea  surface. 

The  received  signal  is  time-averaged  for  a  duration  that  is  long  compared  with 
the  sea  surface  fluctuations.  This  is  a  limiting  case  that  is  not  directly 
reported  in  the  literature;  however  it  is  relatively  easy  to  devise  an 
adequate  expression  to  account  for  this  case.  Only  the  time-averaged  nns 
pressure  is  needed;  there  is  no  need  for  consideration  of  coherent  and 
incoherent  scattering,  or  the  statistics  of  the  time-varying  surface. 

This  report  proposes  simple  equations  for  the  time-average  ms  received 
pressure  from  an  omnidirectional  source  in  the  vicinity  of  the  sea  surface. 


Figure  1  -  Geometry  of  the  Image  interference  effect 
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At  low  frequencies  where  absorption  is  not  significant,  and  neglecting 
refraction  effects,  the  mean-square  pressure  at  a  receiver  in  the  vicinity  of 
a  flat  pressure-release  surface  can  be  written 


2  2  r 

P  r 
ro  o 


L1  L 


2  2 


2yr. 


1  + 


cos  2-rrf  T 
o 


0 

where  r^  and  r 2  are  ray  distances  as  shown  in  Figure  1;  fQ  is  the  signal  fre- 

r2_ri 

quency;  T  *  — - —  ■  the  time  difference  in  arrivals  via  ^  and  r2;  y  is  the 
surface  reflection  loss  coefficient  (can  vary  from  0  to  1).  (In  the  general 
case  where  the  sound  speed  c  is  not  a  constant.,  refraction  will  be  present,  so 
r i  and  r2  will  not  be  straight  lines.)  For  the  sake  of  simplicity,  let  the 
reference  pressure  PQ  *  1  and  r^  *  1. 


Average  over  a  frequency  rand 


It  can  be  shown  (Appendix  A)  that,  for  a  band  of  frequencies  of  bandwidth 
B  and  arithmetic  center  frequency  fn  the  mean-square  received  pressure  may  be 
written  as 


2  yr^  sin  ir$T 
r2  Ttftr 


cos  2itf  T 
o 


(2) 


where  the  factor  may  be  interpreted  as  a  dimensionless  cross-correlation 

TlfT 

coefficient  (with  time  delay  T),  between  the  sound  pressures  from  the  direct  and 
surface-reflected  paths.  This  equation  only  can  apply  when  the  sea  surface  is 
smooth  enough  so  that  the  ocean  wave  height  is  small  compared  to  acoustic  wave¬ 
length,  or  when  the  surface  scattering  is  small.  If  this  is  not  the  case,  a  more 
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general  expression  must  be  obtained  that  takes  account  of  the  sea  roughness 
as  well  as  bandwidth. 


EFFECT  OF  ROUGH  SEA  SURFACE 

The  following  five  points  constitute  a  heuristic  argument  j  arrive  at  a  more 
general  theory  for  the  omnidirectional  transmission  ..oss  as  a  function  of  both 
bandwidth  and  sea  surface  roughness  (wave  height) . 

1.  Most  of  the  discussions  of  surface  loss  (2)  in  the  literature  apply 
to  narrow-beam  sources  and/or  receivers.  Actually  the  value  of  referred  to 
in  this  case  is  a  partitioning  factor  indicating  the  relative  amount  of  sound 
energy  that  is  scattered  out  of  the  specular  direction,  that  is,  the  ratio 
of  specularly-reflected  pressure  to  the  total  incident  pressure.  The  justif¬ 
ication  for  this  interpretation  of  is  that  there  is  very  little  actual  energy 
loss  due  to  reflection  of  sound  from  the  sea  surface.  The  available  mechanisms 
for  such  real  dissipative  surface  loss  are,  in  general: 

a)  Transmission  through  the  surface:  causes  a  very  small  reflection  loss 
of  .005  dB 

b)  Absorption  by  bubbles:  also  estimated  at  about  .005  dB  below  10  kHz 

c)  Absorption  by  organisms:  varies,  but  generally  very  small. 

At  high  frequencies,  (say  50  kHz)  the  loss  due  to  bubbles  may  become  significant 
at  low  grazing  angles  where  the  sound  travels  through  an  extensive  thickness  of 
bubbles.  Therefore,  with  this  exclusion,  we  can  let  1  and  remove  it  from 
equations  1  and  2. 

2.  Since  ship  noise  measurements  are  almost  always  time-averaged  measurements, 
it  is  not  necessary  to  treat  coherent  and  Incoherent  reflected  sound  separately, 

as  is  often  done  In  the  literature  (2).  Equation  (2)  is  expressed  in  a  con¬ 
venient  form  such  that  there  can  be  one  correlation  factor,  g,  which  will  com- 
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pletely  describe  the  effect  of  the  surface  on  reflected  sound,  for  any  degree  of 
surface  roughness. 

3.  Equation  (2)  describes  a  form  of  the  image-interference  effect  as  a 
function  of  bandwidth  B  which  has  the  following  properties: 

a)  For  B«0,  equation  (2)  reduces  to  equation  1,  which  is  the  case  of  image 
interference  from  a  flat  surface. 

b)  For  B»2fQ,  equation  (2)  reduces  to  the  limiting  case  of  wideband  noise: 


which  is  purely  a  function  of  geometry.  This  function  is  shown  in  Figure  2. 

The  ordinate  in  this  figure  is  plotted  in  terms  of  the  increase  in  pressure 

level  (  in  dB)  due  to  the  presence  of  the  surface  path.  The  0-dB  line 

represents  the  free-field  transmission  loss  20  log  rj.  The  abscissa  is 

plotted  in  terms  of  the  dimensionless  geometry  parameter  Va.  (The  levels 

ri 

given  are  strictly  applicable  only  to  the  case  of  isovelocity  water.)  Notice 
however,  that  the  maximum  contribution  from  the  surface  path  is  only  3  dB.  This 
Indicates  why  exact  theoretical  derivations  (such  as  by  integrating  the  results 
for  narrow-beam  theory)  are  unnecessary  and  a  simple  definition  of  correlation 
is  sufficient.  Measurement  precision  and  omnidirectionality  of  real  transducers 
rarely  are  better  than  1  dB. 

4.  The  dimensionless  surface  roughness  parameter  commonly  applied  is  the 
Rayleigh  parameter 

o  2irfoa 

R  -  k  osin  a  •  2tt —  sin  a  - - sin  a 

\> 

where  o  is  the  rms  wave  height,  a  Is  the  grazing  angle  of  the  specular  ray,  and 
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GEOMETRY  PARAMETER 


■P  TEAIT  Xfuino  M013U30 


Figure  2 .  Increase  in  mean  level  of  wideband  noise  signals  due  to 
surface  reflected  sound  as  a  function  of  geometry  (square-law 
detector  assumed) . 
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XQ  is  the  wavelength  of  sound  of  frequency  f^.  In  the  case  of  scattering  from 

a  narrow  beam  of  sound,  a  has  a  unique  value.  This  is  the  case  usually  treated 

in  the  literature.  However,  in  the  case  of  an  omnidirectional  source,  sound  is 

can 

scattered  at  all  angles  from  a  rough  surface.  Nevertheless,  we/retain  the  specu- 
lar-reflectioi^alue  for  a  in  the  latter  case.  The  justification  for  this  is  that 
a  is  always  the  angle  from  which  most  of  the  sound  will  be  reflected  on  a  time- 
average  basis,  and  angles  within,  say,  10  percent  of  a  will  usually  cover  a 
large  portion  of  the  scattered  sound  pressure  (3).  Therefore  the  Rayleigh  par¬ 
ameter  will  be  retained  with  the  above  extension  of  its  definition  implied  in  the 
case  of  an  omnidirectional  source. 

If  the  sea  surface  is  flat,  equation  (2)  applies,  for  any  band  of  frequencies 
centered  at  f^.  But  if  we  assume  a  very  rough  surface,  such  that  R  >  10,  it  is 
clear  that  the  reflected  ray  will  be  uncorrelated  with  reBpect  to  the  direct 
ray,  regardless  of  the  value  of  B.  In  other  words,  even  with  B-0,  equation  (3) 
will  apply  when  R  is  large. 

5.  Both  B  and  R  affect  the  correlation  factor.  It  seems  likely  that  they 
both  affect  it  to  the  same  degree,  because  both  are  dimensionless  ratios  of 
lengths,  the  limiting  values  are  the  same,  and  they  are  both  conservative:  they 
only  affect  the  received  phase;  there  is  no  amplitude  loss.  From  the  above  con¬ 
siderations  we  infer  a  combined  expression  for  the  correlation  factor  of  the 
form  = 


sin  (FBT  +  R) 
TTBT  +  R 


(4) 


It  is  clear  that  the  effects  of  B  and  R  are  additive,  because  when  one  term  is 
small  and  the  other  large,  the  large  one  dominates.  It  also  seems  reasonable  that 
both  B  and  R  have  equal  weight  in  determining  the  correlation  factor,  g. 
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Therefore  a  general  expression  for  the  received  mean  square  sound 

pressure  from  an  omnidirectional  point  source  and  receiver  under  a  rough  sea 


surface  may  be  written 


2r.  sin  (TST+R) 

-  -  co8  2irf  T 

r  -nBT+R  ° 


(5) 


AVERAGE  OVER  A  FINITE-SIZE  SOURCE 

Equation  (5)  may  be  extended  further  to  include  the  effect  of  spatial  aver¬ 
aging  in  a  vertical  dimension  caused  by  a  source  of  finite  size.  This  situation 
occurs,  for  instance,  in  the  case  of  a  cavitating  propeller  on  a  surface  ship, 
where  the  collapsing  cavitation  bubble  is  extended  in  depth  and  radiates  sound  across 
this  extended  region. 

In  this  case  the  depth-averaging  term  has  a  form  that  is  identical  to  the 
Rayleigh  parameter,  except  that  the  rms  wave  height  s  is  replaced  by  the  rms 
vertical  source  width  d: 

Z  -  2  %  d  sin  a  -  2  a  d  fo  sin  a 
X  c 


This  is  the  correlation  factor  for  spatial  averaging  in  the  vertical 
dimension.  Clearly  this  factor  is  additive  with  the  other  factors. 

Therefore,  the  resultant  image  interference  equation  is  obtained  by  adding  all 
three  factors: 


■  sin  (*BT  +  R  +  Z) 

8  “  1TBT  +  R  +  Z 

so  that  the  resulting  image  interference  equation  becomes 
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<?% 


2  2 

P  r 
*o  o 


2 

r2 


2r 


1  sin  (it  BT+R+Z) 
BT+R+Z 


cos  2rr  f  T 
o 


(6) 


This  model  of  the  image  interference  anomaly  has  been  implemented  in  a 
BASIC-language  program  for  the  Hewlett-Packard  9845  computer.  A  listing  is 
included  in  Appendix  B.  The  program  uses  conventional  1/3-octave  frequency 
bands  and  can  accept  up  to  4  receiver  depths  in  a  vertical  array,  with  two 
commonly-used  methods  of  averaging  the  data  from  different  receivers. 
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Appendix  A.  Derivation  of  Image  Interference  Term 
for  a  Band  of  Frequencies 


Equation  (1)  applies  only  to  the  case  of  a  single  frequency  fQ,  whereas  It 
is  more  usual  to  measure  noise  in  a  band  of  frequencies  B  *  f2  -  f In  this 
case  the  third  term  in  equation  (1)  must  be  averaged  over  frequency: 


2tBT 


cos  2nfT  df. 


'fl 


fl+f2 

Defining  fQ  -  — - — ,  0^  «  irf^T  and  ©2  “  irfjT,  we  obtain 


r 

rBT  Is 


2nBT  I  sln  2"£2T  *  sln  2 


2(0 


1 


2m [ 

^7  [’ 


sin  27Tf2T  -  sin  2irf 


"£iT] 

2l] 


sin  202  -  sin  20^ 


2  1 


sin  202  sin  20 2 


] 


] 


■  6^7  [8l“  6  2  co>92  '  ,ln9l  c0,9l] 

1  T  2  2  2  2  i 

-  q  _Q-  [(cos  0j^  +  sin  01)  sin02  cos02  -  (sin  02  +  cos  02>  sinQj^  cose1  J 


.  .  2  q,  cos  0«  -  sin  0_  cos  0,  sin  0,  -  cos“  0-  sin  0-  cos  0„ 

sin  0-  cos  °1  2  v2  V1  7.  1  2 


0^[ - 2 


+  cos  @2  sin  0 


1  sin  02] 


sln  @2  cos  0^  -  cos  02  sin  0 


e2^i 


1  1 
-  ^cos  0^  cos  02  -  sin  0^^  sin  02J 
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COS  (6^+62) 


sin  (02  -0j) 


e2"ei 


s  in  nBT  _ o  rp 

— — — —  cos  2ntT. 
ttBT  o 


12 


18  !  LLOYDS  LLOYD-MIRROR  INTERFERENCE  FOR  1/3-OCTAVE  MERSUREMENTS 

20  !  PAUL  RRVESON  REVISED  JUNE  2,  1983 

30  I  DATA  STORAGE  VERSION  OF  LLOYDP 

40  I  BASED  ON  FORTRAN  PROGRAMS  BICENT  (CCPAP),  PSM  (CCPAM),  AND  OTHER  MORK 
90  !  This  program  generates  anomaly  in  propagation  loss  due  to  surface 

60  I  image  interference  as  well  as  anomalies  due  to  hydrophone  direction- 

70  !  ality,  absorption  loss,  and  detector  type.  Options  are  available 

00  !  to  account  for  averaging  due  to  a  vertical  source  width 

90  i  and  hydrophone  averaging  method  < average  of  dB  levels  or  average  of 

100  !  powers  before  dB  conversion). 

110  !  The  program  does  not  account  for  sometimes  significant  effects 

120  !  due  to  averaging  over  range  during  sample  time,  refractive  anomalies, 

130  !  or  bottom  reflections. 

140  OPTION  BASE  1 

190  PRINTER  IS  16 

160  GCLEAR 
170  PRINT  PAGE 

180  DIM  D2 ( 1 0  > , Freq  <  4 1 ) , P  < 1 8 ) , Asum <  4 1 ) , Dbsum  <  4 1 ) , Db ( 1 8 , 4 1 > 

190  DIM  Fre  <  4 1 ) , A  < 1 0 , 4 1 ) , Ak  yds  <  2 1 ) , X  <  4 1 > , Y ( 4 1 > , P$ ( 6 ) , C$  C  80  3 
200  SHORT  Xsum<34> 

210  DATA  .0093, .0068,. 0089,. 104,. 13,. 164,. 219,. 293,. 419,. 604,. 9, 1.36, 2. 07 
220  DATA  3.14,4.73,7,10.1,14.1,18.8,23.9,29 
230  MAT  READ  Ale  yds 

240  RESTORE  210 

290  Fmax-34  !  UPPER  FREQUEHCY  IS  20,000  Hz 

260  Msuslf*":" 

270  Det-1 

280  PRINT  " PROPAGATION  ANOMALY  DATA" 

290  BEEP 

300  INPUT  "SQUARE  LAU  DETECTORS  ARE  ASSUMED;  IF  LINEAR  ENTER  1",L* 

310  IF  L8»" 1 "  THEN  Det-PI/4 

320  BEEP 

330  INPUT  "ENTER  LENGTH  OF  HYDROPHONE  ELEMENT  IN  INCHES",Hy1 
340  Hyd-Hyl/12 

390  BEEP 

360  INPUT  "ENTER  OCEAN  HAVE  HEIGHT  IN  FEET",Hh 
370  BEEP 

380  INPUT  "ENTER  MEAN  SOURCE  DEPTH  IN  FEET",D1 
390  BEEP 

400  INPUT  "ENTER  SOURCE  VERTICAL  WIDTH  IN  FEET",Dx 
410  IF  Dl>*Dx  THEN  GOTO  460 

420  BEEP 

430  DISP  "ERROR  —  SOURCE  NOT  FULLY  SUBMERGED" 

440  WAIT  1900 

490  GOTO  370 

460  BEEP 

470  INPUT  "ENTER  HORIZONTAL  RANGE  IN  YARDS", Cpay 
400  BEEP 

490  INPUT  "HON  MANY  HYDROPHONES  ARE  AVAILABLE?" , Nhyd 
900  IF  NhydOO  THEN  GOTO  990 

910  BEEP 

920  DISP  "BAD  INPUT,  TRY  AGAIN" 

930  GOTO  400 

940  HAT  B2-<0> 

990  FOR  1*1  TO  Nhyd 

960  BEEP 

970  INPUT  "ENTER  DEPTH  OF  HYDROPHONE  IN  FEET",D2<I> 

900  NEXT  I 

590  I  START  OF  COMPUTATIONS  - 

600  Cpa«3*Cpay 

610  Cpaq*Cpa*2 

620  C-5000  I  Sound  speed  In  feet  per  second 

630  Avgl*l 

640  Avg2«0 

650  DISP  "WAIT  A  MINUTE  —  I'M  THINKING" 

660  FOR  1*1  TO  Nhyd  I - - -  13 - HYDROPHONE  DEPTH  LOOP  - - 


ere  sum>oi+B2<i) 

680  Biff«Bl-B2(I> 

690  Rsq«Cpaq+Biff*2 
700  Rrsq"Cp*q+SuaA2 

710  R>SQR(Rsq)  I  Source  to  hyd.  distance 

720  Rr«SQR<Rrsq>  I  Image  to  hyd.  distance 

730  P(I>«SQR(D1«S2<I>>/'R  !  Geoaetry  parameter  for  surface  reflections 

740  Thetal-PI/2 
750  Theta2-PI/'2 
760  RAB 

770  IF  Cpa>0  THEN  Thetal«ATN(Bi  ffsCpa)  I  Angle  source  to  hyd. 

700  IF  Cpa>0  THEN  Theta2*RTN(Sua^Cpa>  I  Angle  Image  to  hyd. 

790  Ang1e<I)«ABSCThetal*180spi) 

600  H-RsRr 
810  T-(Rr-R)'C 

820  FOR  F-l  TO  Fmax  ! -  FREQUENCY  LOOP - 

830  Freq>10*10M.ie<F-l>> 

840  Fre(F>m0R0UN9CFreq,3) 

850  E1«E2«1  I  Absorption  loss  negligible  below  1  KHz. 

860  Laabd4»C/Freq 

870  IF  F<21  THEN  GOTO  910 

880  A1pha«Akyds<F~20>s60000  i  Divide  by  3000*20  to  give  exp.  loss  per  foot 

890  El»10M-Alpha*R>  I  Absorption  loss  factor  for  R 

900  E2>10A(~R1pha*Rr)  i  Absorption  loss  factor  for  Rr 

910  01 r l-FNDi rec tCThetal , Lambda, Hyd ) 

920  B i r 2-FHB 1 rec  t  <  The t  a2 , Lambda, Hyd  > 

930  Cl»El*81rl 

940  C2«E2*01r2 

950  81na*Sum^R  I  Sine  of  angle  from  image  to  hyd. 

960  B*.2316*Freq  1  Bandwidth  in  hertz 

970  Z1«PI*B*T  I  Bandwidth  Integral 

980  Z2>2*PI*0x*S1na^Laabdal  Integral  over  source  width  (HEH  VERSION) 

990  Z3B2*PI*Mh*Sina>'Laabda!  Integral  due  to  surface  roughness  (Rayleigh) 
1000  X*21+Z2+Z3 

1010  Corr*SIN(X)sX 

1020  Ad ,F>«Cl*2*(M*C2)A2-2*M*Cl*C2*Corr*C08<2*PI*Freq*T)  I  Lloyd  mirror  eqn 
1030  Bbd,F)»10*LGT<Rd,F)*Bet) 

1040  NEXT  F 

1050  Avgl«Avgl*  Angled) 

1 060  Avg2«Avg2*Ang  led) 

1070  NEXT  I 

1080  Avgl«AvglAd<'Nhyd)  i  Average  angle  computed  as  harmonic  mean 

1090  Avg2«Avg2sNhyd  I  Average  angle  computed  as  ordinary  mean 

1100  FOR  F«1  TO  Fmax  !  Compute  hydrophone  data  averages  two  ways 

1110  Asum<F)*0 

1120  8bsun<F)*0 

1138  FOR  I-l  TO  Nhyd 

1140  Asum(F)»Asum<F)*A(I,F)*Bet 

1150  Bbeum<F>»Bbsum<F>*Bbd,F> 

1160  NEXT  1 

1170  Asum<F)>10*LGT(Asue<F>sNhyd)  f  Asum  -  AVERAGE  OF  POHER  VALUES  IN  OB 

1100  Bbsum(F)«Bbsum(F)/'Nhyd  I  Bbsum  -  AVERAGE  OF  BB  VALUES 

1190  NEXT  F 

1200  0I8P 

1210  PlotdonemB 

1220  Store:  1  8T0RE  A8UN  VALUES  ON  BISK  MITH  SUPPORT  BATA 

1230  BEEP 
1240  Yo»no*0 

1250  INPUT  "STORE  COHPUTEB  VALUES  ON  BISK?  <C0NT»N0,  l-YES)", Yesno 
1260  IF  Ye»no«0  THEN  GOTO  Print 
1270  BEEP 

1200  EBIT  "ENTER  HASS  STORAGE  UNIT  SPECIFIER", Hsuslt 
1290  NASS  STORAGE  IS  NsuslS 
1300  BEEP 

1310  INPUT  "ENTER  NEW  FILE  NANE" ,Na»e8 

1320  00m"  "  I  CREATE  SUPPORT  BATA  STRING 
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1338 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1700 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1070 

1000 

1090 

1900 

1910 

1920 

1930 

1940 

1900 

1960 

1970 

1900 

1990 


FIXES  1 
H1*-VAL*(Hyl > 

Hh»«VAL#(Wh> 

Sd#-VAL#(B1> 

Sw««VALt(Bx> 

Hr*-VAL*(Cp«y> 

Hdl*>VALt(B2(l>> 

Hd2*-VAL*(B2(2>) 

Hd3*-VAL*(B2<3)> 

Hd4*-VAL*(82(4>> 

C*-N*»*#l,S*l."HL«"«.H1»tS$t"MH«"l.Wh$l,S**."SB-"tSd*I.S#t"SM«"iSw#68#fHR«"4Hr* 

C*-C»tS*t"HB-"t.Hdl*iS*tHd2*t.S#iHd3*tS$tHd4* 

PRINT  C* 

FOR  K*1  TO  34 
Xsu«(K>*A»um(IO 
NEXT  K 

CREATE  Nwil,l 
ASSIGN  «1  TO  N*»*t 
HAT  PRINT  «i;Xsun 
PRINT  «l;C* 

ASSIGN  tl  TO  • 

PRINT  "LLOYB-HIRROR  SATA  STORES  IN  FILE  "iNwtf 

Prints  !  -  PRINT  OUTPUT  - 

BEEP 

Hc>8 

INPUT  "PRINTEB  SATA  FOLLOWS;  WANT  HARBCOPY?  (l»YES>",Hc 
IF  Hc-1  THEN  PRINTER  IS  0 

PRINT  "PROPAGATION  ANOHALY  VERSUS  NYSROPHONE  BEPTH  ANS  FREQUENCY" 

PRINT 

PRINT  "FROH  FILE  ";Naa«* 

PRINT 

IF  L#< >" 1 "  THEN  PRINT  "SQUARE  LAN  BETECTORS  ASSUHEB" 

IF  L*-"l"  THEN  PRINT  "LINEAR  BETECTORS  ASSUHEB" 

PRINT  "VERTICAL  LENGTH  OF  NYSROPHONE  ELEHENT  •  "}Hyi;"  INCHES" 

PRINT  "OCEAN  NAVE  HEIGHT  -  "JMH;"  FEET" 

PRINT  "HEAN  SOURCE  BEPTH  ■  ";B1J"  FEET" 

PRINT  "SOURCE  VERTICAL  MIBTH  ■  "jBx;"  FEET" 

PRINT  "HORIZONTAL  RANGE  -  ";Cp*yJ"  YARBS" 

PRINT  LIN<3> 

PRINT  "HYBROPHONE  BEPTH  <FEET>  ANB  ANGLE  (BEG. >" 

IHAGE  12X,4<BBBB,4X>, "PHR  AVG.  BB  AVG. " 

PRINT  USING  1730; B2( 1 >lB2(2>yB2(3>fB2(4> 

IHAGE  13Xt6(BB.B,4X> 

PRINT  USING  1750;AngU<l>f Angl •  <2 ) , Angl •  <3) , Angl •  <4 ) , flvgl , Awg2 
PRINT  "FREQUENCY" 

PRINT 

FOR  F»l  TO  10 

IHAGE  DBDBBD . B , 3X , 6 ( SBB . B , 3X > 

PRINT  USING  1800;Fr«<F>t8b<l,F)y8b<2,F>,Bb<3,F>,Bb<4yF>;A*ci»<F>;Bb*u»<F> 
NEXT  F 

FOR  F«ll  TO  Faax 
IHAGE  2X(BBBBBBy3Xy 6(800. B,3X> 

PRINT  USING  1840;Fr*(F>yBb(lyF>yBb(2yF>yBb(3yF>yBb(4yF>;A»w»(F>;Bb»u»(F> 
NEXT  F 

IF  Hc>l  THEN  PRINT  PAGE 

PRINTER  18  16 

Yn«0 

BEEP 

INPUT  "HANT  PLOT?  (l«YE8>",Yn 
IF  YnOl  THEN  GOTO  LMtlln* 

Pit:  I  GENERAL  BATA  ANB  CURVE  PLOTTING,  CARTESIAN  COORB I NATES 

I  GRAPH  SETUP  BATA 
Y«1n«-30 
Ymx«10 
X«1n-l 
XmxMI 
NvtrtH 
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2008  Nhorz*40 

2010  Xlabel*-" FREQUENCY  BAND,  Hz" 

2020  Y 1  zb* 1 4* "ANOMALY ,  dB" 

2030  XspaABS  <  Xb&x-Xm 1 n  > /Nhor z 
2040  Yap- ABS<Yaax-Ya1n> /Avert 
2050  Ygap-ABS<Yaex-Ya1n>e.e05 

2000  Plot graph:  !  GENERAL  GRAPH  SETUP  **************** 

2070  PRINT  PAGE 

2000  PLOTTER  IS  "GRAPHICS" 

2090  GRAPHICS 
2100  FRAME 
2110  SETGU 

2120  LOCATE  12, 118, IS, 93  I  Deflnea  graph  irtt  within  100  X  123  fraae 
2130  SETUU 

2140  SCALE  Xain,Xaax,Yain,Yaax 

2150  LINE  TYPE  3 

2160  GRID  10,5,1,-30 

2170  LINE  TYPE  1 

2100  AXES  l,5,Xain, Ya1n 

2190  AXES  l,S,Xaax,Yaax 

2200  1  LABEL  GRAPH  *#***##*#*#####*##«#♦##♦*** 

2210  DEG 

2220  LDIR  90 

2230  LORG  9 

2240  C8IZE  2.4 

2250  FOR  1-1  TO  41  i  Label  x-axla  nuabera 

2260  HOVE  I, Yaln-Ygap 
2270  LABEL  Fred) 

2200  NEXT  I 
2290  LDIR  0 
2300  LORG  0 

2310  FOR  J-Yain  TO  Yaax  STEP  Yap  I  Label  y-axla  nuabera 

2320  HOVE  Xain.J 

2330  LABEL  J 

2340  NEXT  J 

2350  SETGU 

2360  CSIZE  3.3 

2370  HOVE  65,2 

2300  LORG  5 

2390  LABEL  Xlabel*  I  Label  x-axla  title 

2400  HOVE  3,50 

2410  DEG 

2420  LDIR  90 

2430  LORG  4 

2440  LABEL  Ytabet#  I  Label  y-axia  title 

2450  LORG  1 

2460  LDIR  0 

2470  HOVE  110,1 

2400  C8IZE  2.4 

2490  LABEL  "LLOYD" 

2500  SETUU 
2510  UNCLIP 

2520  Piet data:  I  PL0T8  DATA  FROH  UP  TO  6  FILES  USING  6  DIFFERENT  8YHB0L8 
2530  P*< !>•"♦" 

2540  P#<2>-"*" 

2550  P6<3>«”X" 

2560  P*<4>«"0" 

2570  P0<5>»"0" 

2500  PS<6>-"*" 

2590  FOR  1*1  TO  Nhyd 
2605  SETUU 
2610  LORG  5 
2620  CSIZE  2.4 
2630  FOR  J«1  TO  Faax 
2640  XU>«3 

2650  YO>-»b<l,J> 
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2660  IF  Y<J>>Y»ax  THEN  Y<J>-Yaax 

2670  IF  Y<JXY«ln  THEN  Y<J>-Yain 

2606  IF  X<J>>X»ax  THEN  X<J>-Xaax 

2690  IF  X<  JXXain  THEN  X<J>-Xain 

2709  HOVE  X<J>,Y<J> 

2710  LABEL  P*(I> 

2720  NEXT  J 
2730  8ETGU 

2740  HOVE  40410*1,90 
2750  LABEL  P#<I> 

2760  HOVE  40410*1,95 
2770  LABEL  D2U> 

2780  NEXT  I 
2790  Finish:  I 

2S00  MAIT  4000 
2010  BEEP 
2020  Copy-0 

2830  INPUT  "NANT  HARDCOPY?  < 1 -YES Copy 

2840  IF  Copy-1  THEN  DUHP  GRAPHICS 

2850  P1otdon«-l 

2060  PRINT  PAGE 

2870  Last  1  in*:  I 

2880  END 

2090  I 

2900  DEF  FND1r*ct(Th*t4,Laabda,Hyd> 

2910  !  NYDROPHONE  VERTICAL  DIRECTIVITY  FUNCTION 

2920  Ar g- ABS  <  P I *Hyd*S I N  <  Th* t  a>  sLaabda) 

2930  Diract-1 

2940  IF  Arg>0  THEN  Dir*ct-SIN<Arg>sArg 
2950  RETURN  Diract 
2960  FNEND 
2970  END 
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